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We consider lattice gas automata where the lack of semi-detailed balance results from 
node occupation redistribution ruled by distant configurations; such models with non- 
local interactions are interesting because they exhibit non-ideal gas properties and can 
undergo phase transitions. For this class of automata, mean-field theory provides a cor- 
rect evaluation of properties such as compressibility and viscosity (away from the phase 
transition), despite the fact that no //-theorem strictly holds. We introduce the no- 
tion of locality - necessary to define quantities accessible to measurements - by treating 
the coupling between nonlocal bits as a perturbation. Then if we define operationally 
"local" states of the automaton - whether the system is in a homogeneous or in an in- 
homogeneous state - we can compute an estimator of the entropy and measure the local 
channel occupation correlations. These considerations are applied to a simple model 
with nonlocal interactions. 

Keywords: Lattice gas automata, //-theorem, entropy, correlations. 

1. Introduction 

Various lattice gas automaton (LGA) models lacking semi-detailed balance 
(SDB) have been proposed because they exhibit interesting properties for in- 
vestigating phenomena which are not accessible to LGAs satisfying SDB. LGAs 
violating SDB locally have been constructed for their operational value as Possi- 
ble candidates for the simulation of high Reynolds number hydrodynamicsEratJ. 
Such systems were subsequently reconsidered as a paradigmatic formulation 
of driven systemsu. The implications and consequences of the lack of semi- 
detailed balance in lattice gas automata are non-trivial, in particular for what 
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concerns their stationary distribution function and the very existence of an 
equilibrium distribution for such systems which exhibit thermodynamic prop- 
erties and transport properties correctly described by mean-field theory. This 
somewhat paradoxical situation has been the subject of consideration in the 
recent literatureou where the question was raised whether the asymptotic state 
of systems without SDB can be qualified as an equilibrium state. Certain LGAs 
without detailed balance were shown to approach a stable uniform state referred 
to as a non-Gibbsian equilibriumu with nonlocal spatial correlations, in contrast 
to the classical Gibbsian state of thermal equilibrium. 

The problem has not been examined in depth for the case of LGAs where the 
lack of SDB results from node occupation redistribution ruled by distant config- 
urations; such models with nonlocal interactions (NLIs) are interesting because 
they exhibit non-ideal gas properties^ and can undergo phase transitional For 
this class of LGAs, mean-field theory is also seen tixprovide a correct evaluation 
of properties such as compressibility and viscositytffl (away from the phase tran- 
sition) and static correlation functions are in accordance with those observed in 
real fluidsH. It has even been suggested that a generalized version of the SDB 
condition could be applicable to recast LGAs with NLIs such that a generalized 
-H-theorem would holdO. The situation calls for clarification, starting at the 
level of the basic statistical mechanical description of LGAs. 

2. Statistical Mechanics 

We consider a set S of Boolean variables (bits), S — {s(^,i = 1, . . . , bV} with 
S(i) G {0, 1} (here bV is simply an integer; b and V will be specified subse- 
quently), and an update rule: T\ : S — > 5, with A a set of parameters, allowing 
for example T to be drawn probabilistically among a set of rules. This update 
rule is set according to certain constraints which are chosen on the basis of phys- 
ical requirements. S and T define the lattice gas automaton. The constraints 
on T are of two types: 

• Geometric constraints lead to a spatial representation of the LGA, and im- 
pose certain symmetries on this representation. This essentially amounts 
to defining a set of relations {(i, <p(i)) ,i = 1, • ■ • , bV} between indices of 
elements of S. With these relations, the update rule T\ may be factor- 
ized into several steps, one of which is simply the copying of the value of 
bit i onto bit tp(i); this step is called propagation and the corresponding 
operator a is denoted by V . An important feature of propagation is that 
it is a mere correspondence between indices of elements of S, indepen- 
dently of the values taken by these elements. To each <p(i) we associate 
a geometric vector ct{\ in a suitable space, and these vectors (which may 
include the null vector) are grouped into b equivalence classes. The vectors 

a This operator may itself be stochastic^, but here we restrict ourselves to the deterministic 
case. 
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representative of these classes will be denoted by Cj (with i = 1, . . . , b); 
they must satisfy symmetry constraints in order to obtain a consistent 
geometric representation. 

• Conservation laws provide physical content. For example, denote by 
7i(j) the value taken by the variable s^y, then define a quantity N[S] = 
Si=i n (i)i which is identified as the number of particles, and demand that 
N{T\S] = N[S]. A similar conservation is generally required for linear 
momentum p = n (i) c i- 

We factorize the transformation T\ as: 

T x =VoU x (1) 

where V is the propagation as defined above. Now, through the operator U\, 
it may happen that one can define locality . If U\ is such that there exists a 
decomposition of S into subsets s.( x ) m such a way that the restriction of U\ 
over s.( x ), denoted by U\\ s . {x) , is completely defined and is an endomorphism: 

U\\s. {x) ■ -» > V(x) , 

— i.e. if the values taken by the bits of s.( x ) after application of U\ are entirely 
defined by their values before this application — then the set of indices that 
define the subset s.^ x ) is called a node] the values taken by the bits of s.m 
define the state of the node. Given the above definition of b, V is now seen to be 
the number of nodes of the automaton universe and b the number of channels 
per node. Then one may construct a representation of the bits that belong to 
a given s.( x ) as 'particles' and 'holes' spatially located on the spot defined by 
x. The action of U is then local and U is called a collision operator (and will 
be denoted by C). We call L-subsets the subsets s.( x ) thus defined. Note that 
it is not mandatory that a LGA have such a property; for example the models 
defined in section 3 have an update operator that cannot be decomposed strictly 
into propagation and collision operators. 

2.1. A Liouville H -theorem 

We denote by T the set of all possible sets S. We associate to each universe- 
state S a probability P(S) in a Gibbs ensemble and define a transition prob- 
ability Ass' between universe-states, such that the system obeys a Chapman- 
Kolmogorov equation: 

P{S') = Y,A ss >P{S),{S'er). (2) 

ser 

The following hypothesis on A, the semi-detailed balance (SDB) condition: 



-Ass> = 1, 

ser 



(3) 
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suffices to prove an iJ-theoremE2l; then the global entropy 

H = -J2P(S)laP(S) (4) 
ser 

does not decrease under the action of U. Furthermore, as noted above, the 
propagation operator produces a deterministic redistribution of the bits without 
introducing or removing information; consequently V does not modify the value 
of if. 

For all practical implementations, |T| is finite, but so large that no mean- 
ingful sampling is realizable; as a result the quantities P(S) and Ass> are not 
accessible to significant measurements. 

2.2. Locality and nonlocality 

If the updating operator T is such that it admits a decomposition into a propa- 
gation operator V and a local collision operator C, the local subsets st x \ define 
subspaces 7( x ) of T, with V = ®^ l(x)- If we assume that the 7( a: )'s are isomor- 
phic, then there exist local states s such that the description of the automaton 
in terms of the s's is equivalent to its description in terms of S. Those states 
(of the nodes) are independent of each other for the local operator, so that we 
can define a local probability p{st x \) in a (smaller) Gibbs ensemble. If the local 
dynamics is described with a transition probability matrix_/l ss ' , then under the 
condition X^e 7 A ss ' = £ 7), a local H -theorem holdalJ and the quantity 

H x ) = - ^2p[ s (x) = s]lnp[s (a .) = s) = -^p(s (a .))lnp(s (a .)) (5) 

does not decrease under the action of C. 

In this local description, the action of V is unimportant as it just provides 
at each iteration a fresh configuration of the bits of sr x \ with probability p(s) 
by definition; this probability is such that h does not decrease. The action of V 
is of importance when one considers the connection between the local entropy 
(||) and the global entropy (^) . If the propagation does not produce correlations 
between the nodes, then the global entropy H is entirely determined by the 
local entropy h: 

fl- = $>(*). (6) 

X 

The existence of a notion of locality provides us with objects (the S( x )'s) 
which are accessible to measurements; then we can define the local number of 
particles n(x) and the local distribution function p(x) by averaging n(x) over 
the appropriate Gibbs ensemble: p(x) = {n(x)}. We may choose the s's as fun- 
damental dynamical objects to which "physical" properties are then associated. 
In usual LGAs, this is accomplished by imposing constraints such as "mass" 
and "momentum" conservation. It is only because the update rule allows a 
definition of locality — i.e. the universe is a set of independent subsets — that 
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such constraints have physical justification since isolated objects can now be 
considered. The local constraints are chosen according to the problem to be 
investigated and provide physical significance to the LGA model. 

An important point is that the decomposition may be further pursued: we 
can describe the dynamics in terms of local space- velocity occupation rii{x) - 
i.e. the value taken by the bit S( x )j, with x = 1, V, and i = 1, b - and of 
local single-particle space- velocity distribution functions p[s^i = 1] = fi(x) = 
/(x, Cj). At this level of description an entropy is defined by: 

b 

h,(x) = -Y^\fi(x)]nMx) + fi(x)]n Jtfx)] , (7) 
»=i 

(where fi(x) = 1 — fi(x)) which, under the Boltzmann ansatz, is simply equal to 
h(x) (H). The second term on the r.h.s of @ stems from the correlations between 
"particles" (1-bits) and "holes" (0-bits) because of the Boolean nature of the 
variables. Omitting this term amounts to neglecting the important contribution 
of the "particle - hole" correlations to the entropy. 

If the operator U does not decompose 5 into independent, disjoint sub- 
sets, then we have no consistent definition of locality. In particular, we cannot 
meaningfully define states s and probabilities A ss i , and we have no definition 
of objects to which Statistical Mechanics can be applied: the semi-detailed bal- 
ance condition is void, and no "local" ii-theorem can exist. An example is the 
class of automata with "nonlocal interactions" (see section 3). We will however 
argue that some models have update rules that allow for a weaker definition 
of locality, i.e. there may exist a collection of disjoint subsets of S, the bits of 
which are strongly coupled through 14, but only weakly coupled to the bits of 
other subsets. We will call these subsets £-subsets. Then we can operationally 
define these subsets as objects of Statistical Mechanics, even if they do not 
contain the full dynamics of the automaton. This is equivalent to ignoring all 
dependencies and correlations between "nodes" . If the coupling between "non- 
local" bits can be treated as a perturbation, then we may be able to extract 
significant information on the dynamics via "local" quantities. However some 
basic elements necessary to establish the validity of mean-field theory (mainly 
the iJ-theorem) are then absent, and care must be taken in defining and using 
"standard" quantities. 

3. A Simple Model with "Nonlocal Interactions" 

An example of a system where semi-detailed balance does not exist is the LGA 
with "nonlocal interactions" introduced by Appert and ZaleskiQ, subsequently 
analyzed by Gerits el all, and generalized by Tribel and BoonEl. Here we con- 
sider an utmost simplified version of the automaton which nevertheless exhibits 
all essential features and physical properties of the original model. 

The LGA is composed of a set of bits (typically 393216) and an update rule T 
as described in Section 2. This rule admits a decomposition into two operators; 
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one of them is the propagation operator V which copies the value contained in 
S(j) onto the bit labeled The actual computation of ip(i) is performed by a 
computer routine which is described elscwheroO. The computation procedure 
exhibits the following features: it defines the topology of the space used for the 
representation of the state of the automaton on the surface of a torus and it 
yields six classes of equivalence of the (i, <p(i))'s to which we associate six vectors 
Cj; for symmetry reasons these vectors are chosen to be coplanar and normalized. 
The spatial representation of the lattice therefore has a toroidal topology and a 
planar geometry, the geometry of a triangular lattice with hexagonal symmetry. 

The second operator U which together with V, composes the full updating 
operator is itself decomposed into two parts: 

U=U 2 oUi. 

U2 groups the bits of S by subsets of six, and re-shuffles the bits within the 
subsets. The grouping is such that, for any i in a given ^-subset s/ x \, ip(i) 
is in another ^-subset s/ x y, for the re-shuffling, we use the rules of the FITP- 
I automaton (the collision outcome is governed by a Boolean parameter )tL3. 
We denote by sr Xt j\ the j-th bit of the subset sr x y, since there is a one-to-one 
correspondence between the indices (i) and the indices (x,j), we may denote 
by ip(x,j) the index of the bit where V copies the bit (x,j), and by (p k (x,j) 
the index where (x,j) should be copied by applying k times the operator V . 
Finally, we will be interested in the ^-subset to which ip k (x,j) belongs; its index 
will be x( x i j> 

The operator U\ re-shuffles the bits that belong to different ^-subsets, ac- 
cording to the following procedure: 

• draw a random number D, with equal probabilities among {1,2,3} (cor- 
responding to the three axial directions of the lattice) ; 

• for each subset sr x \, draw a random number r according to a given proba- 
bility distribution (which may be degenerate); then, write the value n( x ,j) 
taken by the bit sr x j\; 

• if ; s (x,D+l) = 1; s (x,D+2) = 0, S( x ( x ,£»,r),0+1) = 0; an d S( x (x,£),r),£»+2) = T 

then exchange these values such that becomes 1 and vice-versa; 

• repeat with D + 4 and D + 5; 

• repeat the last two steps with D replaced by D + 3. 

All the values of D + E(E = 1,...,5) are mapped onto (1,...,6). The set of 
parameters A in U\ is formed by the Boolean parameter of the collision, the 
direction D, and the distance r (the set of distances). 

This model exhibits a phase transition when r exceeds a given value (or an 
average value computed over a given probability distribution); away from the 
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Figure 1: Estimated entropy, Eq.(8), as a function of time (t > t* , see text). 
FHP automaton with strictly local collision rules and SDB (solid line), and 
LGA with NLI (interaction distance = 5 lattice units): (i) entropy estimated 
from the channel distributions fi (indistinguishable from full line), and (ii) en- 
tropy estimated from the configuration probability p(s,t) (dashed line). All 
cases: average density / = 0.1, lattice size 256 x 256 nodes, t* = 500; time 
unit=automaton time step. 

transition, the macroscopic properties of the LGA are characterized by well- 
defined coefficienta3'El. 

3.1. Entropy 

We may operationally define "local" states s(x,t) (^-subsets) of the automaton 
(at time t), and, if the LGA is in a homogeneous state, we can evaluate the 
occurrence probability p(s, t) of a configuration s by measuring the occurrence 
frequency u(s, t) over the whole lattice at each time step. A second measure of 
p(s,t), closer to its definition, is the occurrence frequency at a given "location" 
of the automaton, over a large number of realizations; at each run the automa- 
ton is initialized independently with a given set of macroscopic constraints. This 
(time-consuming) procedure is necessary if the automaton is not in a homoge- 
neous state at all times. The occurrence probability changes with time, and so 
does the quantity: 

h cst .{t) = -J2"(^t)^v(s,t)- (8) 

s 

In an automaton with strictly local rules and satisfying SDB this quantity 
will increase monotonously towards a maximum value; in a nonlocal automaton, 
this increase is not guaranteed. 
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Figure 2: Example of transition probabilities yielding a seemingly valid SDB. 
Left boxes: initial configurations (arrows indicate particles). Right boxes: out- 
put configurations (after non-local interaction and collision) with the corre- 
sponding probabilities, assuming that the particles on the "pairing nodes" are 
randomly distributed (p is then simply one-third of the average density). Con- 
sidering the sum of the transition probabilities with an output configuration 
where the particle is in the upper right channel (shown as dot-dashed arrow), 
the sum is clearly equal to one. 

Consider a FHP type LGA subject to a constraint realized by imposing 
a systematically oriented output configuration for binary collisions; then after 
some time f*, the constraint is relaxed by restoring the usual rule with output 
configurations equally distributed along the three directions. During the first 
phase the entropy decreases from its initial value, then at time t* it starts 
increasing and levels off at its maximum value. This example is illustrated in 
Fig. |l| for LGAs with and without NLIs. For a LGA with strictly local rules and 
satisfying SDB, we observe monotonic increase of h est as expected; for the LGA 
with NLIs, the entropy estimated using the configuration probability p(s) shows 
non-monotonic increase and stabilizes at some value below the plateau value 
obtained for the LGA without NLIs. Alternative interpretations are possible: 
(i) if h is a valid quantity and is correctly measured — i.e. v is a good estimator 
of p — then we conclude that SDB is violated and the i?-theorem does not 
hold; (ii) the local states are here ill-defined, and we conclude that h is not 
a meaningful quantity and that the //-theorem at this level of description is 
irrelevant. Now LGAs without SDB can also be viewed as operational models 
to simulate non-equilibrium systemsE^I; but LGAs with NLIs would be driven 
systems without either boundaries or any external "field" , which renders the 
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Figure 3: The entropy estimator, evaluated using the occurrence frequency 
v(s,i) of a "local" configuration at a given node of the lattice, for 150 time- 
steps and 10000 realizations. The system is initialized with a step-function in the 
density field which is then relaxed; the estimator is measured on a node located 
at a position along the initial step. Parameters for the simulation: distribution 
of interaction distances between 1 and 8 lattice units, with p(r) oc r _1 ; initial 
density: / = 0.1 on one half of the lattice, / = on the other half; lattice size 
256 x 256 nodes. No phase transition occurs for these values of the parameters. 

"non-equilibrium constraints" rather unphysical. 
3.2. Transition probability 

Using the weak definition of locality (section 2.2), we define as above "local" 
states of the automaton. Then a transition frequency F ss * can also be defined 
and we can measure F ss > either by averaging over the whole "lattice" (if the 
state is homogeneous) or by averaging over many realizations (Gibbs ensemble). 
Again we must stress that this is not a strictly correct estimator of the transition 
probability A ss > which does not exist in automata with nonlocal rules where the 
transition frequency depends on the configuration of the whole set S (i.e. on 
the density field over the whole lattice); the transition frequency measured in 
an inhomogeneous state differs from its value obtained in a homogeneous state. 

Consider a LGA model with "nonlocal" interactions where we can measure 
the transition frequency F ss i . Here we are only interested in "local" configura- 
tions, but we can make no a priori assumption either on the direction of the 
interaction or on the state of the "pairing" nodes. We then define a "mean- 
field" transition probability A™/', evaluated by considering an arbitrary node 
and assuming that the configurations of all the other nodes are given by a prob- 
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Figure 4: Nonlocal channel pair occupation probabilities, evaluated on nodes 
separated by interaction distance r in LGA with NLI as described in the text. 
The probabilities Pij(r, cd) are measured along a chosen direction D (here D — 
0) when the interaction between channels labeled i and j occurs along that 
direction. The black dots are the superposition of data for i — ... 5. The 
density field is homogeneous; parameters of the simulation: / = 0.1, r = 5, 
lattice size 256 x 256 nodes. 

ability distribution p(s). h If we assume that the bits are all uncorrelated and 
that each channel has occupation probability / (the average density per chan- 
nel), any configuration has probability p(s) = / JV ( S )(1 — f)l b ~ N (. s )] > where N(s) 
is the number of occupied channels in configuration s. Then the interactions are 
equally probable along any direction, and since the interaction rules have the 
same symmetries as the underlying lattice, each momentum change can occur 
with the same probability as the reverse change. Therefore in this mean-field 
picture, the collision rules are such that SDB is satisfied, and so is the complete 
update rule, i.e. 

£A?/=1. (9) 

Note that @ does not mean that the matrix elements A ss t are identical in 
systems with and without nonlocal interaction; only the sums over all initial 
states are equal to one in both cases; for an example, see Fig. g. However the 
time behavior of the entropy estimator is different in the two types of system; 
when NLIs are present the entropy estimator is not a monotonously growing 

F ss i is a statistical measure of the unknown quantity A gs i, while A m 'f' is a theoretical 
evaluation of this quantity under the mean- field hypothesis. 
c Hcre even detailed balance is satisfied. 
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Figure 5: Local channel correlations (Cy(0) = (nj(x)nj(x)) — f 2 ). Parameters 
of the simulation: / = 0.1; lattice size 256 x 256 nodes; interaction distance 
r = 5, 10; for comparison the open circles refer to the LGA without NLI. 

function, see Fig.[l] (homogeneous case) and Fig.|| (inhomogeneous case). Direct 
measurement of the transition frequencies reveals indeed that the mean-field 
approximation is very poor (we find significant deviations to Eq.(^)) in accor- 
dance with the fact that the iJ-theorem does not hold. We conclude that the 
mean-field transition probability is not an appropriate quantity to correctly de- 
scribe the dynamics of LGAs with nonlocal rules. Furthermore the reasoning 
does not include the interaction distance, so omitting one of the most essential 
features of the model. 

3.3. Boltzmann ansatz 

In LGA theory, the Boltzmann ansatz allows to establish the explicit connection 
between different levels of description of the _ff-theorem (hence the existence of 
a known equilibrium), and provides a way to obtain the explicit expression 
of the channel occupation distributions /, (x, t) . The first point is essentially 
of theoretical importance; the second point has crucial consequences for the 
computation of the properties of LGAs. 

Nonlocal interactions occur only between nodes with certain configurations, 
producing output configurations deterministically. If we denote the direction 
of interaction at time t by D, then only particles on the pairs of channels 
{D + 1,D + 2} and {D + 4,13 + 5} (modulo 6) can interact, thereby ex- 
changing their momentum. As a result, a configuration where pairs of chan- 
nels {{D + i, D + i ± 1}) are occupied by "converging" particles becomes more 
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probable than a configuration with "diverging" particles. Since the direction 
of interaction is randomly chosen at each time step, this effect is not visible 
in single-node configurations, but must be expected to show up in nonlocal 
channel-pair correlations CV,- (r, cry ) defined by 

Pij(r,c D ) = (nj(c)nj(x + rc D )) 

= / 2 +a,(r,c D ), (fO) 

with / the average density. The Boltzmann ansatz assumes that 

C i , j (r,c D ) = f(l-f)6 ij 6 r0 . (11) 

This is indeed the case in a LGA with SDB; however, in automata with non- 
local interactions, or more generally without SDB, correlations arise. For the 
present model, one expects to observe correlations between channels oriented 
along a 60-deg. angle, since only particles on channels with neighboring in- 
dices do interact. Figure |J shows the correlations measured in a simulation; 
the results are clearly in agreement with the predicted effect. The LGA Enskog 
formalismQ can be generalized to incorporate higher orders in the correlations; 
a comparative analysis between theory and simulation results will be presented 
elsewhere. 

Nonlocal interactions are sources of correlations. Once they are created, 
these correlations propagate, and, at sufficiently low density, they may survive 
until the involved particles reach the same node (following their initially con- 
verging trajectories). The results given in Fig. ^ show indeed that nonlocal 
correlations (induced by NLIs) create local correlations (i.e. between adjacent 
channels on the same node). The existence of such correlations explains the 
difference between the values of the entropy estimated via the single-particle 
distribution functions fi, and of the entropy estimated via the configuration 
probability p(s) (see Fig. [j]): local correlations create an entropy deficit. Indeed 
it is a general property that the entropy computed from the full statistics of a 
set of Boolean variables is lower than, or equal to the entropy computed from 
the set of individual statisticsE3: correlations contain information, and therefore 
reduce the entropy. 

4. Comments 

We have presented a description of lattice gas automata using a rigorous defini- 
tion of locality. The basic objects to be used in LGA statistical mechanics and 
in simulation measurements are defined in terms of the operator which updates 
the bits of the automaton. The analysis introduces a distinction between (i) 
LGAs where a strict notion of locality exists, in which case the semi-detailed 
balance condition is well defined (whether satisfied or not), and (ii) LGAs with 
a weak notion or no notion of locality, where the states of the node do not de- 
termine the full dynamics, and where semi-detailed balance loses meaning. For 
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the first class of models the question of the existence of a local equilibrium is 
quite relevant. In automata of the second class, the existence of a global equi- 
librium is, in general, not related to the existence of a local equilibrium, and 
the Boltzmann hypothesis is invalid at all ranges. We have proposed a weaker 
definition of locality and of local states which can be used in appropriate cases; 
with this lax definition, mean-field analysis can be justified for the computation 
of thermodynamic and transport properties of the lattice gas. 
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